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EXPERIMENTAL  AND  COMPUTER-SIMULATION 
STUDY  OF  HIGH-VELOCITY  IMPACT 
DEFORMATION  IN  PREHEATED  MATERIALS 


ABSTRACT 


Ws  have  developed  a new  “reverse  gun”  technique  for  studying  high-strain-rate 
deformations  in  preheated  materials.  In  this  procedure,  a rigid  plate  is  fired  at  a stationary 
heated  cylindrical  rod,  in  contrast,  the  standard  rifle  procedure  fires  an  unheated  rod  into  a 
stationary  plate.  The  advantage  of  our  reverse  gun  method  is  that  the  rod  can  be  preheated 
and  the  temperature  monitored.  We  validated  our  technique  by  comparing  both  methods  in 
experiments  conducted  at  room  temperature,  and  we  obtained  new  deformation  data  for 
aluminum  and  copper  rods  heated  to  725  K.  We  observed  good  agreement  between  the 
measured  shape  of  the  recovered  rods  and  the  shapes  generated  by  hydrodynamic  code 
calculations  using  a modified  Steinberg-Guinan  constitutive  model  and  a Griineisen  equa- 
tion of  state._ 


INTRODUCTION 

The  U.S.  Air  force  has  expressed  interest  in  an  air-delivered  penetrator  that  can  pierce  tank  armor. 
One  design  developed  for  this  type  of  weapon  incorporates  a metal-lined  high  explosive  lens  that,  when 
detonated,  ejects  a slug  of  metal  (“self-forging  fragment”)  that  forms  a rod-shaped  penetrator  during  the  First 
moments  of  its  flight.  The  penetrator  retains  this  shape  until-  it  strikes  the  target. 

Hydrodynamic  computer  codes  used  to  design  such  a device  and  calculate  its  performance  require 
input  information  concerning  the  behavior  of  materials  at  elevated  temperatures  and  at  very  high  strain  rates 
(<  > 10^  see*').  The  dependence  of  flow  stress  on  temperature  and  strain  rate  exhibited  by  potential  penetrator 
materials  was  therefore  of  particular  interest,  and  we  initiated  a feasibility  study  at  LLL.  to  generate  these  ex- 
perimental data.  At  the  same  time,  we  began  a program  of  computer  simulation  in  order  to  interpret  the  ex- 
perimental results.  We  describe  here  these  preliminary  experiments  and  our  computer  simulations. 


MATERIALS 

The  new  experimental  technique  was  initially  tested  on  rods  made  from  6061-T6  aluminum  and 
oxygen-free,  high-conductivity  (OFHC)  copper  obtained  from  laboratory  stock.  These  rods  received  no 
special  heat  treatment  prior  to  these  tests. 

The  materials  of  primary  interest  were  two  types  of  electrolytic-tough-pitch  (ETP)  copper  furnished 
by  Eglin  Air  Force  Base.  Axially  and  transversely  oriented  rods  were  obtained  from  billets  designated  1 14  mm 
(Type  A)  and  152  mm  (Type  B)  ETP  copper.  Metallurgical  details  and  mechanical  properties  are  reported 
elsewhere.*  The  billets  were  machined  into  cylindrical  rods  6 mm  diam  by  30  mm  long  that  were  annealed  at 
575  K for  one  hour  and  cooled  in  air. 
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EXPERIMENTAL  METHOD 

Two  methods  were  used  to  study  high-strain-rate  plastic  deformation.  Data  were  obtained  for  sam- 
ple rods  at  room  temperature  with  a 7.62  mm  rifle  in  the  manner  described  by  Taylor^  and  by  Wilkins  and 
Guinan^(Fig.  I).  In  this  system,  a moving  cylindrical  rod  with  flat  ends  strikes  a rigid  wall.  Assuming  that  the 
rod  decelerates  as  a rigid  body  and  that  plastic  deformation  occurs  at  the  impact  surface,  we  can  equate  the 
rate  of  decrease  of  rod  length  L with  the  velocity  U: 

^ - u (1) 

The  deceleration  of  the  rod  is 

d^L  _ dU  , Yq 

— — TT r ’ 

dt2  dt  PqL 


where  Yq  is  the  flow  stress  (assuming  an  elastic-perfectiy  plastic  model)  and  pgis  the  initial  density  of  the  rod. 
Combining  Eqs.  (I)  and  (2)  and  integrating  yields 


" ■ 28n(Lf/Lo)  ’ 

where  Lg  and  Lf  are  the  initial  and  Hnal  lengths  of  the  rod.  The  flow  stress  can  thus  be  estimated  through 
measurements  of  density,  projectile  velocity,  and  the  initial  and  final  lengths  of  the  rod.  However,  because  the 
assumptions  that  plastic  flow  occurs  only  at  the  impact  surface  and  that  flow  stress  is  constant  are  only  ap- 
proximations, an  accurate  relationship  between  the  flow  stress  and  the  final  length  of  the  projectile  is  much 
more  complex  than  Eq.  (3).  The  flow  stress  is  generally  a function  of  several  variables  and  can  be  obtained 
only  by  comparison  of  experimental  data  with  hydrodynamic  code  calculations. 

If  the  velocities  of  the  colliding  bodies  are  interchanged  so  that  a projectile  equivalent  to  a rigid  body 
strikes  an  unrestrained  rod,  Eqs.  (l)-(3)  are  unchanged.  This  second  method,  called  the  “reverse  gun",  is  il- 
lustrated in  Pig.  2.  The  advantage  of  this  system  is  that  the  stationary  rod  may  be  heated  prior  to  impact  which 
allows  the  relation  between  plastic  deformation  and  temperature  to  be  studied. 


FIG.  2.  Prqectile  and  target  configuration  for  “reverse  gun”  experiments. 


The  reverse  gun  experiments  were  conducted  using  a com  pressed -gas  gun  shown  schematically  in 
Fig.  3.  The  launch  tube  was  2.44  m long  with  a bore  63.5  mm  in  diameter.  The  impact  parts  of  the  projectile 
consisted  of  a disk  of  alumina  (AI2O3)  25.0  mm  in  diameter  and  12.0  mm  thick  with  density  3.73  Mg/m^.  The 
alumina  was  backed  by  an  aluminum  disk  25.0  mm  in  diameter  and  15.0  mm  thick.  These  parts  were  bonded 
to  the  front  face  of  a 63.5-mm-diam  piston  as  shown  in  Fig.  2.  Projectile  velocity  was  determined  using  time- 
interval  counters  that  measured  the  time  between  pulses  generated  by  the  short-circuiting  of  “whisker  pins” 
separated  by  calibrated  distances  along  the  bore. 

The  cylindrical  targets  were  held  as  shown  in  Fig.  4.  A complete  momentum  exchange  was  assured 
by  designing  the  massive  projectile  and  target  array  so  that  complete  deformation  of  the  target  rod  occurred 
before  any  significant  interaction  between  the  projectile  plate  and  the  peripheral  target-holding  apparatus 
could  take  place;  the  needle-pointed  support  screws  were  notched  so  that  they  fractured  easily.  Variation  of 
the  plate  thickness  showed  that  target  deformation  was  independent  of  this  parameter. 

The  rods  were  heated  to  about  725  K using  infrared  radiation  supplied  by  a furnace  located  as  shown 
in  Fig.  3.  The  furnace  was  designed  so  that  it  could  be  dropped  from  the  projectile  path  a moment  before  the 
gun  was  fired.  Rod  temperature  was  monitored  at  the  instant  of  impact  by  means  of  a digital  readout  from  a 
thermocouple  cemented  to  the  heated  rod. 

The  momentum  trap,  which  allowed  soft  recovery  of  the  target  rods,  was  a lOO-mm-diam  by  700- 
mm-long  cylinder  of  polyurethane  (p  at  2 Mg/m^)  with  a 15-mm-diam  by  200-mm-long  hole  drilled  in  the 
front  end.  Following  impact  and  full  target  deformation,  the  relatively  unrestrained  target  rod  was  propelled 
through  the  1 5-mm-diam  hole.  This  configuration  segregated  the  rod  from  the  target-holding  and  projectile 
debris.  The  polyurethane  at  the  bottom  of  the  hole  decelerated  the  target  rod  without  deformation. 

The  validity  of  the  reverse  gun  data  for  6061-T6  aluminum  at  room  temperature  was  established  by 
comparison  with  new  7.62  mm  rifle  shots  and  with  rifle  shot  results  from  Ref.  3 (Fig.  5). 

Some  additional  room  temperature  tests  were  performed  on  rods  of  OFHC  copper  fabricated  from  a 
single  billet  of  LLL  stock.  The  reverse  gun  data  were  compared  with  data  from  6.0-mm-diam  by  30-mm-long 
cylinders  of  OFHC  Cu  that  were  placed  in  plastic  sabots  and  fired  in  a 7.62  mm  rifle  at  an  alumina  disk 
12.5  mm  thick  backed  by  32  mm  of  steel.  Figure  6 compares  the  results.  Disagreement  with  data  from  Ref.  3 is 
believed  to  result  from  actual  material  differences,  probably  preshot  work  hardening. 
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FIG.  4.  Target-rod  holding  ap- 
paratus. The  rod  is  supported  by 
needle-pointed  screws.  The  entire  sup- 
port system  was  spring-loaded  to  allow 
for  thermal  expansion. 


FIG.  S.  Comparison  of  reverse-gun 
data  (triangles)  with  7.62  mm  rifle 
data  performed  by  us  (circles)  and 
from  Ref.  3 (smooth  curve)  for  room- 
temperature  6061-T6  aluminum. 


FIG.  6.  Comparison  of  reverse-gun 
data  (triangles)  with  7.62  mm  rifle 
data  performed  by  us  (circles)  and 
from  Ref.  3 (smooth  curve)  for  room- 
temperature  OFHC  copper.  The  dis- 
crepancy between  our  data  and  those  of 
ReC  3 represent  a significant  dif- 
ference in  material  preparation. 


0.1  0.2  0.3 

U — km/s 


RESULTS 


The  experimental  data  are  summarized  in  Tables  1-4.  The  cross-check  data  for  606I-T6AI  and 
OFHC  Cu  are  displayed  in  Figs.  S and  6.  The  plots  for  isotherms  at  295  K and  725  K for  ETP  copper  are 
shown  in  Fig.  7 and  constitute  the  primary  experimental  results;  the  strong  temperature  dependence  of  the 
length  ratio  for  ETP  copper  is  clearly  evident.  The  295  K isotherm  also  demonstrates  that,  as  far  as  these  tests 
are  concerned,  there  is  no  significant  difference  in  Lf/Lo  between  the  axially  and  transversely  oriented  rods,  or 
between  the  type  A and  type  B materials. 

These  results  clearly  demonstrate  the  feasibility  and  validity  of  the  reverse  gun  technique. 


6 


TABLE  1.  Data  summaiy  for  ETP  copper  (T  = 295  K). 


Shot 

number 

TypeCu* 

Initial  den.’-V 
Mg/m^ 

Projectie  velocity 
(U),  km/f 

Pins  Optics 

Final  length/ 
Initial  length 
(Lj/Lo) 

23 

A-a 

8.907 

0.148 

b 

0.806 

24 

A-a 

8.921 

ai81 

b 

0.675 

25 

A-a 

8.928 

0.316 

0.307 

0A18 

26 

A-a 

8.919 

0.277 

0.258 

0.492 

27 

A-a 

8.935 

0J57 

0.340 

0.332 

28 

A-a 

8.928 

0.367 

0J54 

0.325 

29 

A-t 

8.918 

0.311 

0J20 

0.373 

30 

B-a 

8.916 

0.332 

0.321 

0.383 

31 

B-a 

8.918 

0.225 

0.213 

0.597 

32 

B-a 

8.911 

0.109 

b 

0.812 

33 

B-a 

8.904 

0.126 

b 

0.777 

34 

B-a 

8.902 

b 

0.394 

0.280 

35 

B-a 

8.911 

0.349 

0.352 

0.333 

*A  = 114-fflin  baiet,  B = lS2-min  bllet,  a « rod  cut  axially  from  bilet, 
t = rod  cut  tianfveiaety  from  bilet. 

^’No  lecoid. 


TABLE  2.  Data  summaiy  for  preheated  ETP  copper. 


Shot 

number 

Type  Cu* 

Initial  density 
(p^,  Mg/m^ 

Temperature 

m.  K 

Projectie 
velocity 
(U),  km/s 

Final  length/ 
Initial  length 
(L,/Lo) 

Comments 

52 

A-a 

8.930 

718 

0.188 

0.555 

53 

A-a 

8.928 

765 

0.237 

0J91 

Large  peripheral 
cracks 

55 

A-a 

8.928 

703 

0.145 

0A90 

56 

A-a 

8.929 

721 

0.210 

0J07 

Bent  off-axis 
during  recovery 

57 

A-a 

8.927 

724 

0.206 

0.529 

Sight  petiphoral 
ctacU^ 

58 

A-t 

8.913 

733 

0.19 

0.574 

U obtained  from 
kndiog  curve 

60 

A-t 

8.929 

717 

0.142 

• 0.718 

61 

A-t 

8.947 

727 

0.211 

0.516 

62 

A-t 

8.935 

721 

0.115 

0.769 

63 

A-t 

8.937 

718 

0.0487 

0.993 

64 

A-t 

8.929 

727 

0.191 

0J62 

*A  > 1 14  mm  bilet,  a ^ tod  cut  axially  from  bUet,  t = tod  cut  tranfveiaely  from  bUet. 
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TABLE  3.  Data  summary  for  OFHC  copper  (T  = 295  K). 


Shot 

Initiil  densiW 

Projectile 

Final  length/ 

Cor^menU 

numbei 

(Pq),  Mg/m^ 

velocity 

Initial  length 

(U),  km/5 

(W 

38 

8.930 

0.177 

0.784 

Reverse  gun  data 

43 

8.907 

0.194 

0.764 

Reverse  gun  data 

45 

8.912 

0.233 

0.688 

Reverse  gun  data 

48 

8.912 

0.0895 

0.933 

7.62  mm  rifle  data 

SO 

8.907 

0.189 

0.774 

7.62  mm  rifle  data 

51 

8.930 

0.276 

0.5% 

7.62  mm  rifle  data 

TABLE  4.  Data  summary  for  6061-T6  aluminum. 


Shot 

Initial  density 

Initial  temperature 

Projectile 

Final  length/ 

Comments 

number 

(Pq),  Mg/m^ 

(T),  K 

velocity 

Initial  length 

(U),  km/s 

(W 

18 

2.685 

295 

0.252 

0.848 

7.62  mm  rifle  data 

19 

2.691 

295 

0.257 

0.841 

7.62  mm  rifle  data 

20 

2.687 

295 

0.361 

0.717 

7j62  mm  rifle  data 

21 

2.651 

295 

0.406 

0.663 

7.62  mm  rifle  date 

39 

2.708 

295 

0.20 

0.903 

Reverse  gun  data 

40 

2.708 

295 

0.358 

0.735 

Reverse  gun  data 

41 

2.710 

635 

0.194 

0.770 

Reverse  gun  data 

42 

2.707 

655 

0.354 

0.397 

Reverse  gun  data 

COMPUTER  CALCULATIONS 


To  conduct  hydrodynamic  computer  code  simulations,  we  need  a complete  constitutive  mode!  that 
incorporates  a pressure  equation  of  state  (EOS)  and  expressions  for  the  shear  modulus  G and  (low  stress  Y. 

Since  the  material  in  our  experiments  is  not  heated  above  the  melting  point,  we  are  interested  in  an 
EOS  for  the  solid  phase  only.  Specifically,  we  choose  the  Mie-Griineisen  model, 

P(V,T)  - P(V,0)  = 2^  [E(V,T)  - E(V,0)]  , (4) 


where  P is  the  pressure,  V is  the  volume,  E is  the  internal  energy,  and  7 is  the  Gruneisen  parameter.  Since  the 
zero  Kelvin  curve  P(V,0)  has  not  been  measured,  we  must  compute  it  from  another  reference  curve,  namely, 
the  Hugoniot.  Thus 


Ph(V)  - P(V.O)  = [E„(V)  - E(V.O)) 

where 

Eh(V)  = ^ Ph(V)  (Vq  - V) 


r 


(5) 


» 
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A T = 295  K,  Type  A - a 

0 T = 295  K,  Type  A - 1 

□ T = 295  K.  Type  B - a 

O T a 725  K,  Type  A - a 

• T a 725  K,  Type  A - 1 


U — km/s 

FIG.  7.  Experimental  ETP  copper  iaotherms  for  295  K and  approximately  725  K showing  linear  behavior  over 
the  meainred  velocity  range.  Tbm  curves  indicate  that,  for  the  same  impact  velocity,  a heated  rod  is  ^gniflcantly 
dMrter  after  iaspact  than  a rod  at  room  temperature.  The  greater  deformation  is  attributed  to  decreased  flow 
stress  at  the  hi^r  temperature. 


I 


Assuming  a linear  shock  velocity-particle  velocity  relation 


Us  = C + SUp  . 


we  obtain 


i 

r 


p„(V)  = -£sSh— 
(1  - Sx)2 


where 


X = I - — . 


A simple  and  moderately  accurate  Gr'tineisen  parameter  is 


r(V) 


I'oV 


Equation  (5)  now  becomes  a diiTerential  equation  for  Eo(V)  > E(V,0).  A power  series  expansion  in  x is 
suitable  for  Eg,  since  only  small  compressions  (less  than  10%)  are  expected  in  our  experiments.  Thus 

EqOO  = Ego  + Egi  x + Eg2  x2  + Egj  X^  + Eg4  X^  , 

P(V.O)  * Pg(V)  = - 3Eg/aV  . 

Assuming  a constant  heat  capacity  3R  for  the  solid,  we  derive  the  following  coefficients: 

Egg  = - 900R  , 

Egi  “ToEoo  • 

®02  “ 2 * “^0  ®0c)  • 

Egj  - ^(4SC2  +73  Egg)  . 

Eo4  “ ^ (-27gSC2  + 18S3c2  + Egg)  . 

For  hydrodynamic  code  calculations,  we  use  the  independent  variables  m - (p/pohl  tuid  e « E/>g 
We  need  expressions  for  the  pressure  P(/t,e)  and  the  temperature  T(M,e).  For  the  compressions  found  in  the  ex- 
periments, it  was  sufficient  to  take  the  pressure  expansion  to  and  the  temperature  expansion  to  /t*.  Thus 
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POi.  e)  = Po(m)  + [e  - CoOi)] 


= PoC2m  + PoC2  |^(l  -I|)  + 2(S  - 1)]  m2 

+ PoO  [2^1  - “lljcs  - 1)  + 3(S  - 1)2]m2  + To®  • (6) 

and 

T(M,e)  = [^e/po  - Eqo  ' EqI#'  " (-^01  + ^02^^  ~ (^01  " ^^02  * Eo3)<^^ 

- (-Eq,  + 3Eo2  - 3Eo3  + Eo4)m^]  /3R  . (7) 

We  now  consider  models  of  the  shear  modulus  (G)  and  flow  stress  (Y).  First,  we  need  an  expression  for  the 
melting  temperature  as  a function  of  volume,  because  G and  Y will  drop  abruptly  to  zero  at  the  melting  point. 
We  use  the  Lindemann  melting  law,^  which  works  well  for  metals  like  A1  and  Cu: 


TmW  = T„o  exp  (27oX)  (1  - x)2/3  • 

The  flow  stress  Y is  a scalar  quantity  and  must  be  compared  with  a scalar  invariant  of  the  stress  tensor. 
Speciflcally,  the  von  Mises  yield  criterion  is  used; 


(oj  - 02)2  + (02  - 03)2  + (03  - a,)2  * 2Y2  , 

where  the  values  of  a are  the  principal  stresses.  If  the  stress  invariant  is  less  than  2 y2,  the  motion  is  elastic;  if  it 
is  greater,  then  the  motion  is  plastic. 


For  G and  Y,  we  use  the  Steinbeig>Guinan  (SG)  model^ 
fof  T < T„  , 

G * Co  + b ♦ h (T  - 300)  j . 

Y « Yo  [ 1 ♦ <l(cP  + ^f)]"  [1  ♦ b ♦ hCT  - 300)j  , 

where 

Yo[l  sY.„; 

fo,  T > T.  , 

G - 0 , 

Y - 0 . (9) 


II 


In  this  case. 


V - — > 

Po 

gP  = equivalent  plastic  strain  = f i>  * , 

■ ■*!  [(;?  - '5)'  * (^5  - ^5)’  * US  - . 

and  ( P is  the  initial  equivalent  plastic  strain.  The  constants  b,  h,  /9,  n,  Gq,  and  Y,„ax  ‘trc  tabulated  for  each 
material;  the  values  of  Yq  in  Eqs.  (3)  and  (8)  arc  not  the  same,  however.  No  explicit  strain-rate  dependence  has 
been  assumed  for  Y in  these  calculations. 

The  SG  model  predicts  that  the  shear  modulus  increases  with  pressure  and  decreases  with  tem- 
perature, dropping  discontinuously  to  zero  at  T > T,„.  The  flow  stress  model  uses  the  same  coefficients  b and 
h for  pressure  and  temperature  dependence  since  the  behavior  of  Y is  expected  to  be  very  similar  to  that  of  G. 
In  addition,  a work-hardening  function  dependent  on  the  coefficients  0 and  n allows  Y to  increase  smoothly 
with  plastic  strain.  The  effect  of  work  hardening  is  given  an  upper  limit  Y„,ax- 

For  simulation  of  impact,  we  used  HEMP,  a two-dimensional  Lagrangian  hydrodynamic  code.^  At 
the  end  of  each  time  step,  the  state  of  strain  and  the  energy  in  each  zone  were  calculated.  These  data  were  used 
to  compute  T and  P and  then  Y and  G for  use  in  the  subsequent  step. 

Our  objective  was  to  compare  the  hydrodynamic  calculations  with  experimental  data  over  a range  of 
impact  velocities  and  initial  temperatures  and  to  determine  how  well  the  SG  constitutive  model  reproduces  the 
final  length  and  shape  of  the  rods.  We  considered  for  this  comparison  6061 -T6  aluminum  and  ETP  copper;  a 
complete  listing  of  coefficients  for  these  two  metals  is  given  in  Table  S.  These  calculations  are  expensive  in 
terms  of  computer  time,  and  therefore  we  ran  only  a few  comparisons  with  experiments. 

Each  calculation  was  set  up  initially  as  an  axisymmetric  rod  with  initial  values  of  p,  e,  and  velocity  U. 
The  front  face  of  the  rod  was  in  contact  with  a rigid  boundary  that  was  perpendicular  to  the  rod  axis.  The  ini- 
tial configuration  with  7 X 26  zones  is  shown  in  Fig.  8.  As  the  rod  struck  the  boundary,  the  equivalent  stress 


TABLE  5.  Constitutive  model  and  EOS 
panmeten  for  ahunhuun  and  copper. 
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(6061-T6) 
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1.99 

b 
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2.0 
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0.32*' 

ai2* 

ai23^ 
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FIG.  8.  Initial  conflgurntion  of  the 
rod  in  HEMP  calculations.  The  rod  b 
a 6.4-nim-diani  by  3l.77>nini-long 
cylinder  divided  into  7 radial  and  26 
axial  xones.  The  rigid  wall  is  at  the  left 
boundary. 


exceeded  Y in  the  front  zones  and  these  zones  deformed  plastically  in  the  radial  direction.  As  stress  waves 
moved  back  and  forth  along  the  rod,  the  forward  momentum  of  the  rod  diminished  gradually,  and  a wave  of 
plastic  deformation  moved  away  from  the  impact  surface.  At  the  same  time,  heating  and  work  hardening  oc- 
curred. Plastic  deformation  eventually  ceased  and  the  rod  underwent  elastic  oscillations  of  small  amplitude. 
This  process  for  Al  and  Cu  was  completed  in  40  to  100  Msec,  depending  on  the  initial  conditions. 

Our  approach  to  Fitting  the  computer  simulation  to  the  experimental  data  evolved  through  a series  of 
calculations.  We  considered  Yq  and  h to  be  adjustable  parameters.  Yg  was  determined  by  comparing  simula- 
tions with  experiments  at  low  velocity  on  room  temperature  rods  in  the  spirit  of  Eqs.  ( l)-(3).  Under  these  con- 
ditions, the  effects  of  heating  due  to  impact  are  negligible  and  the  term  that  contains  h is  therefore  unimpor- 
tant. 

The  curve  of  Lf/Lg  vs  U for  aluminum  at  room  temperature  was  approximated  reasonably  well  by 
the  SG  model  if  Yg=  0.32  G Pa  (rather  than  the  0.29  GPa  value  reported  in  Ref.  5).  For  the  shots  with 
preheated  rods,  we  used  the  EOS  to  determine  the  values  of  the  independent  variables  m and  e corresponding 
to  P = 0 and  the  desired  value  of  T.  It  is  signifleant  that  the  temperature  coefficient  h = -6.2  X 10~*,  given  in 
Ref.  5,  is  too  small  to  simulate  the  preheated  shots  and  predicts  a value  of  Lf  that  is  too  large.  Several  trials 
demonstrated  that  a better  value  is  h ==  -I.S  X 10~^.  Using  these  new  parameter  values,  our  calculation  fits  the 
experimental  room  temperature  curve  adequately,  and  it  also  fits  shot  41,  with  T = 615  K,  rather  accurately. 
An  attempt  to  fit  shot  42  failed  because  the  numerous  rezones  required  by  the  extreme  flattening  of  the  front 
zones  impaired  the  accuracy  of  the  calculation. 

T o this  point,  we  have  only  discussed  our  attempts  to  predict  the  final  length  of  the  rod.  A sensitive 
test  of  the  constitutive  model,  and  especially  of  the  effect  of  work-hardening  on  flow  stress,  is  the  shape  of  the 
rod  after  impact.  Figure  9 compares  the  shape  of  the  three  aluminum  simulations  and  experiments. 

A series  of  HEMP  runs  in  which  U » 0.189  kiih/s  and  Yg  was  varied  showed  that  the  room  tem- 
perature copper  data  of  all  types  was  covered  by  the  range  0.12  < Yg  < 0.25  GPa.  This  extent  of  this  range  of 
values  appears  to  reflect  the  different  work  hardening  experienced  by  the  various  types  of  copper  prior  to  our 
experiments. 

We  were  primarily  interested  in  fitting  the  room  temperature  and  preheated  ETP  data  and  found 
that  a value  Yg*  0.125  GPa  provided  an  adequate  fit  of  the  room  temperature  ETP  curve.  An  attempt  to  fit 
shot  52,  which  was  prejheated  to  718  K and  had  an  impact  velocity  of  0.188  km/s,  once  again  showed  that  the 
SG  temperature  coefficient  was  too  small.  Changing  the  value  of  h from  -3.8  X 10^  to  -9.2  X 10"^  yielded 
adequate  agreonent  with  the  measured  final  length. 

Two  additional  long  runs  (to  100  ms)  were  performed  to  check  the  calculations  against  shots  26 
(room  temperature)  and  61  (preheated).  Fracture  at  the  rim  of  the  impact  surface  in  shot  26  prevented  ac- 
curate comparison  of  the  measured  shape  with  the  calculations  for  this  surface.  Comparisons  of  measured  and 
calculated  shapes  for  shots  26,  52,  and  61  are  shown  in  Fig.  10. 
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FIG.  9.  ComptrisoB  of  HEMP 
calcolatioBS  (sBMOtk  curves)  with  ex- 
pcriBMBtal  BieasureuMuts  (circles)  of 
the  profiles  of  deformed  alumhium  rods 
from  shots  39,  40,  sad  41.  The  R 
(radial)  aad  L (axial)  axes  are  aot  to 
scale,  (a)  U *=  0.20  km/s,  T 295  K; 
(b)  U - 0.358  km/s,  T = 295  K;  (c) 
U - 0.194  km/s,  T > 635  K. 


The  agreement  between  experiment  and  calculation  it  good,  considering  that  only  two  of  the 
parameters  in  the  SG  model  were  varied.  The  parameters  in  the  work-hardening  term  as  given  in  Ref.  S thus 
appear  to  be  well  chosen.  Work-hardening  effects  arc  responsible  for  the  inflection  observed  in  the  profllet  of 
the  deformed  rods.  Because  the  experiments  appear  to  have  a consistently  steeper  profile  in  the  inflection 
region  than  the  calculations,  a better  fit  of  the  calculated  curves  to  experiment  might  be  obtained  by  adjusting 
the  work-hardening  parameters  n and  0 slightly. 

The  SG  model  assumes  that  Yq  is  a fixed  property  of  the  material  and  that  any  previous  work  hard- 
ening is  included  in  the  initial  plastic  strain  rl|.  We  have  instead  set  fi|  ■ 0 and  have  assumed  that  materiab 


FIG.  10.  Comparisoa  of  HEMP 
calculatkias  (smooth  carves)  with  ex- 
pcriascatal  measuremcots  (circles)  of 
the  profiies  of  deformed  ETP  copper 
rods  from  shots  26,  52,  aad  61.  The  R 
(radial)  aod  L (axial)  axes  are  aot  to 
scale,  (a)  U > 0.277  km/s, 
T = 295  K;  (h)  U - 0.188  km/s, 
T>7I8K:  (c)  U>=  0.211  km/s, 
T - 727  K. 
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with  different  work-hardening  histories  have  different  Yq  values.  In  either  case,  Yg  or  fij  are  treated  as  ad- 
justable parameters. 

In  these  calculations  we  changed  h,  the  coefTicient  of  temperature  dependence  in  the  shear  modulus 
and  flow  stress,  in  order  to  fit  the  final  length  of  the  preheated  rods.  This  procedure  was  not  quite  correct, 
since  the  temperature  dependence  of  the  shear  modulus  had  been  measured  and  it  should  have  remained  at  the 
measured  value.  In  future  calculations,  the  temperature  coefficients  of  G and  Y should  be  varied  indepen- 
dently. 
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CONCLUSIONS 

The  reverse  gun  technique  effectively  measures  the  temperature  dependence  of  material  properties  at 
high  strain  rates.  In  principle,  the  material  can  be  heated  up  to  its  melting  point  and  measurements  taken  over 
the  entire  temperature  range  in  which  the  material  is  solid.  This  measurement  capability  should  prove  useful  in 
fundamental  studies  of  solid  mechanics  as  well  as  in  hydrodynamic  code  simulation  of  self-forging  fragments. 

The  significant  theoretical  conclusion  of  this  work  is  that  the  SG  constitutive  model,  with  some 
minor  modifications,  can  furnish  an  adequate  representation  of  plastic  deformations  in  aluminum  and  copper 
over  a range  of  impact  velocities  and  initial  temperatures.  | 

Further  experiments,  in  which  a greater  range  of  materials,  impact  velocities,  and  temperatures  are  l 

studied,  would  allow  a more  exact  check  on  the  dependability  of  the  SG  model.  At  high  impact  velocities  and 
high  temperatures,  new  phenomena  such  as  fracture  and  shock  melting  could'be  studied  and  models  could  be 
developed  for  their  interpretation. 
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